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I. INTRODUCTION 



Ultracold atomic quantum gases have been established 
as a clean and tunable test ground for many-body physics 
They allow to mimic condensed matter, but offer 
also opportunities to study aspects of many-body physics 
that are hard to address in other systems. An impor- 
tant example for the latter is the broad and widely un- 
explored subject of many-body dynamics. In this pa- 
per we investigate theoretically the spinor dynamics of 
a Fermi gas far away from equilibrium. We consider a 
one-dimensional system of repulsively interacting spin- 
1/2 fermions confined in a harmonic trap. The initial 
state is created out of the spin-polarized equilibrium by 
rotating the spins spatially into a spiral configuration 
(as previously done for a Bose condensate and pro- 
posed for strongly interacting fermions for the purpose 
of probing the Stoner transition Q). We show that al- 
ready weak interaction, like in the experiment described 
in Refs. [1,01: is sufficient to induce a robust spin seg- 
regation. The effect builds up on times that are long 
compared to the oscillatory motion of the atoms in the 
trap. It can be explained as a consequence of an effec- 
tive, dynamically created long-range interaction that we 
obtain by integrating out the rapid oscillatory dynamics 
in the trap. Within the framework of a semiclassical the- 
ory, the effective interaction is isotropic in phase space. 
The fact that away from equilibrium already weak inter- 
action can cause a noticeable spin segregation contrasts 
with the equilibrium physics of the system. For example, 
the spin segregation of itinerant ferromagnetism, possi- 
ble signatures of which have recently been observed in 
a cold-atom system Q, requires strong interparticle re- 



ulsion as well as higher dimensions. The experiment 
^ has inspired also theoretical work on the dynamics of 
strongly-interacting spin-1/2 fermions, e.g. Refs. 0, 
In this paper we stick, however, to the regime of weak in- 
teraction, where the fermionic cold-atom system does not 
suffer from effects originating from the coupling to molec- 
ular two-body bound states such as dissipative particle 
losses [1] or non- universal scattering properties @ . 

In the following two sections we will first introduce the 
system and describe the semiclassical mean-field theory 
that we use to simulate its dynamics. In section IIVI we 
then present our numerical results, predicting a dynami- 
cal spin segregation. By integrating out rapid oscillations 
in the trap, in section|V]we derive an effective description 
for the dynamics that explains this finding in an intu- 
itive fashion as a consequence of a dynamically induced 
long-range interaction. Before we close with conclusions, 
experimental signatures are discussed. 



II. SYSTEM 

We consider a gas of fermionic atoms of mass M having 
two relevant internal states, m = =ti i- The gas is 

not necessarily quantum degenerate but sufficiently cold 
and dilute to interact via low-energy s-wave scattering 
only. Consequently, the interaction between two particles 
is captured by a pseudo potential g' Pm{mi,m'^m2^{'''i ~ 
r2) with (mi, 7712) and {m'i,m'2) denoting the spin state 
before and after scattering, respectively. Here 
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with unity matrix em'm and vector of spin-1/2 ma- 
trices Sm'm, projects onto the antisymmetric spin sin- 
glet state two scattering particles have (due to Fermi 
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statistics and the symmetric s-wave state). The term 
1 



4 ^m^mi ^m'^m^ 



corresponding to the so-caUed exchange interaction, 
gives rise to spin-spin coupling. The coupling con- 
stant g' = Anh^as/M is proportional to the singlet s- 
wave scattering length as, characterizing the actual in- 
teratomic potential. With this, we can write down 
the Hamiltonaian H' = /dr V;^, (r)/ij„,,„(r)?/'m(r) + 

|- V'mi (^)V'm;, {'r)Pm{m^,m'^m2'4'nl2 (r")V'mi (r). Re- 

peated spin indices imply summation, '0,„(r) is a 
fermionic field operator, and h'^,j^{r) — — ^jjV^em'm + 
V^,^{r) denotes the single-particle Hamiltonian contain- 
ing the potential V^,^{r) = V{r)e,n',n + B{r)-Srn',n that 
can be decomposed into a spin- independent term V'{r) 
and an effective magnetic field B{r) acting on the spin. 

We are interested in the regime where the dynamics 
is reduced to one spatial dimension and consider a har- 
monic confining potential V'{r) = V{x) + V±{y,z) = 
^M[uj^x'^ + Ld^iy^ + z^)] with the a tight transversal 
confinement Ltj± being large compared to other relevant 
energy scales such as the initial temperature or chemi- 
cal potential. Therefore, the particles basically occupy 
the transversal single-particle ground state. Moreover, 
we assume that the magnetic field varies in x-direction 
only, B{r) ~ BqCz + B{x). A possibly present constant 
part BqBz of the magnetic field will be dropped in the 
following, i.e. we are working in a spin frame rotating 
around the z-axis. Introducing a dimensionless descrip- 
tion in units of the longitudinal trap [energies, lengths, 
momenta, and times are given from now on in multiples 
of hjj, {MLo/hY^/"^, {Mfuvf-^^, andw^^, respectively] we 
arrive at the Hamiltonian 



(2) 



+ 1 / da;'(/'^,^(a;)'0^^(a;)P,„'^™i^„,^™2?/>™2(a;)Vimi(a;) 



for the one-dimensional problem. Here g ~ 2hio^as x 



and 



hin' ra{x) c^^x^'ni' va ~\~ ^m' m{x^ , (3) 



(4) 



By swapping field operators and indices {ip^ , -ip^ , 



—tb^ ,tjj^, —ih^ iip^ ,) the interaction can be sim- 

plified to Hint = 9 Jdxip^{x)ip^{x)ipj^{x)ipi{x) refiecting 
that by Pauli exclusion only fermions of opposite spin in- 
teract. The resulting spin coupling - parallel spins avoid 
repulsion - is expressed more clearly, however, in Eqs. ([T]) 
and (12). 



III. SEMICLASSICAL DESCRIPTION 
A. Equations of motion 

We study the system's dynamics in terms of the single- 
particle density matrix 



nm'm{x\x) = tr{p?/i^,(a;')V'm(x)} 



(5) 



with density operator p. It evolves in time according to 
i/5 = [-ff,/5], giving 



yn,-n,'m{x' ,X) = {['li)'^^,{x')ilrn{x),H]). 



(6) 



using cyclic permutation under the trace. While for 
non-interacting particles the r.h.s. of this equation reads 
Kik{x)nra'k{x' ,x) - hkm'{x)nkmix',x), the interaction 
Hint leads to quartic expectation values that we decom- 
pose like 



in order to get a closed equation for n,„',„(x', x). By 
Wick's theorem this decomposition is exact for the initial 
state considered here, being an equilibrium state of a 
quadratic Hamiltonian modified only by the spiral spin 
rotation generated by another quadratic Hamiltonian. At 
later times it corresponds to the time-dependent Hartree- 
Fock approximation that is suitable for weak interaction 
and leads to the non-linear equation of motion 

i'T-m'rn(2; ,2;) = hJ^f,{x)nm' k{x , x) — /l^/ (x )nkm{x , x) . 

(8) 

The Hartree-Fock Hamiltonian h™^,^{x) = h,nim{x) + 
V™,^j^{x) comprises the mean- field potential V^,\„{x) = 
Vni{{x)em'm + Bni{{x) ■ s^'m where 

Vmi{x) ^gno{x), 
B^i{x) = -2gnix), (9) 
X, x) and spin 



with particle density no{x) — e 
density n{x) = Sm'mnm'm{x, x). 

In a next step, we introduce the Wigner function 



iix,p) 



1 

2^ 



d^e 



,{x-i/2,x + £,/2). (10) 



Using Eq. this phase-space representation of the 
single-particle density matrix can be shown to evolve like 

(e.g. m 

1 1 ' ^ 

'W,n'7n{x,p) = -pda:Wm'mix,p) + T X] ^ (^'^^'^p) 

Q = 

'x[Vmk{y)Wm'k{x,p) - {~)°'Vkm'iy)wkmix,p)] (11) 

y=x 

with V,n'm{x) = V,n'm{x) + V^,^^{x) . The relation 

(x, x) ^ dp {x,p) (12) 
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connecting the spatial densities entering the mean-field 
potential to the Wigner function closes this equation of 
motion. We employ a semiclassical approximation to the 
motional degrees of freedom by truncating the infinite 
series after a = 1. This is justified as long as the poten- 
tial Vm'm{x) varies slowly compared to the single-particle 
wave lengths involved. It is, thus, particularly suitable 
for sufficiently hot or dense systems, with cither the ther- 
mal or the Fermi wave length small. Moreover, the trun- 
cation is exact for harmonic or linear potentials Vm'm{x). 
It gives 

Wyn'yn = -pdxWm'm + {dxV)dpWm'm 
-i-B • {SmkWm'k ~ Skm'Wkm) 

+ -{dxB) ■ dp{SmkWm'k + Skm'Wkm) (13) 

having the form of a Boltzmann equation, lacking the 
collision integral and augmented by a coherent spin- 
dynamics, cf., e.g., [TH and references therein. On the 
r.h.s., the four terms describe diffusion, spin- independent 
acceleration, coherent spin rotation, and spin-dependent 
acceleration, respectively. It is instructive to introduce 
the real- valued density and spin Wigner functions 



wi){x,p) = e 
wix^p) = s 



m' m^ra' m 
m' m^m' m 



{x,p). 



(14) 



Their time evolution is determined by (e.g. Ref. [Tl[) 



Wo 



w 



( - pdx + xdp + {dxV^i)dp)wQ 

-t- {dxB + dxB^f) ■ dpW, 

{-pdx + xdp + {dxVmf)dp + {B + JB„if) x )w 



+ ^{dxB + dxBn,f)dpWo 



(15) 



where B ^ only during the preparation while during 
the time evolution to be simulated B = 0. 

The quantum and fcrmionic nature of the system en- 
ters into the equations of motion ()15p in different ways: 
through the initial state {wq,w), via the coherent spin 
dynamics, and with the structure of the mean-field inter- 
action ([9]) stemming from the projection on spin-singlct 
scattering ((T|) for Fermi-statistics. Equations ([T5|) de- 
scribe the coUisionless regime of weak interaction (im- 
plicitely assumed when introducing the mean-field ap- 
proximation [HI). The coUisionless regime is opposed to 
the hydrodynamic regime where collisions constantly re- 
store a local equilibrium of the momentum distribution 
such that the state is described by density and velocity 
fields depending on x only. Both regimes can be found 
in the very same system: The one-dimensional collision- 
less description by the set of equations ([15]) is also valid 
when the transversal confinement uj±_ is not tight enough 
to freeze out transversal motion completely (as assumed 
here), but still tight enough to ensure quick equilibration 
along the transversal directions instead . The results 
presented in this paper are, therefore, also valid in this 



quasi-one-dimcnsional regime. The regime in between 
the coUisionless and hydrodynamic limit is captured by 
adding a collision integral to Eqs. (|15p tending to restore 
thermal equilibrium on a scale Tcor. Considering the ef- 
fect of collisions along the longitudinal direction x be- 
comes necessary when considering stronger interactions 
and longer times scales as we do. For such a regime, it 
has been proposed to observe the spin-wave instability 
predicted by Castaing in a quantum gas p^ . 



B. Initial off-equilibrium state 

Initially, the system is prepared in its spin-polarized 
equilibrium state, with the spins pointing in x-direction, 
temperature T, and chemical potential ^. One has 



^o(a:,p)-— {cxp[-(-p2 + -x2-/i)] +l} (16) 

and w = (1,0, 0)*wo/2 (with), or Wm'm = wo/2. The 
zero-temperature chemical potential, the Fermi energy 
for just one spin-state, simply reads Ep = N with to- 
tal particle number N = Jdxno{x) = JdxJdpwo(x,p). 
This description of the initial state is accurate, since the 
spin-polarized gas is non-interacting and the semiclassi- 
cal approximation exact for a harmonic trap. In a next 
step, at time t = 0, during a short preparatory pulse a z- 
polarized magnetic field gradient is applied, captured by 
B(x) = qxS(t)ez. A spin spiral of wave length As = 2Tr/q 
is created, while 'Wo{x,p) is still given by Eq. (|T6)) one has 



w{x,p) — ( cos(q.T), sin(ga;), O) 



t woix,p) 



(17) 



01 Wrn'm 

{x,p) = exp[iqx{m - m')]wo{x,p) /2. With a 
simple spin rotation, we have prepared a state far from 
thermal equilibrium. Apart from (i) having created a 
rather artificial spiral spin configuration P7)) (favorable 
neither with respect to energy nor entropy), we have also 
increased the number of available single-particle states 
from one spin state to two. The latter has two conse- 
quences: (ii) the phase-space density configuration (jl6p 
is far from being thermal (for half the number of particles 
per spin state having the same kinetic energy as before, 
a thermal distribution is characterized by a lower chemi- 
cal potential and a higher temperature), and (iii) we have 
suddenly introduced interaction to the system. The com- 
bination of (i) and (iii) will lead to robust dynamical spin 
segregation. 



C. Semiclassical versus mean- field dynamics 

When integrating the time evolution for many 
fermions, the semiclassical phase-space equations (|15p 
are usually much easier to treat numerically than the 
Hartrce-Fock mean- field equations. ([5|), even though the 
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interaction is non-local in the former, 

B^f{x) ^ -2g Jdpw{x,p). (18) 

This is exemplified by om- state ([TS]) and (|17p : In phase 
space, it has a linear extent A ~ ma,x[fi, T) [we de- 
fine A = 2max(V2iV,2Vr)], while it varies on the 
scale 6 ~ min(T/A, A^) stemming either from the for- 
mer equilibrium [roughly estimating ^|^gf ~ ^A with 
e = \{x^ + P^) and s = x,p\ or from the induced spin 
spiral. For a reasonable phase-space representation one, 
thus, requires a grid of just more than (A/(S)^ ~ N 
points. On the other hand, the real-space single-particle 
density matrix njn'm{x' , x) varies on the much shorter 
length (5' '--^ A~^ in each argument, calling for more than 
> grid points. 

IV. SIMULATION OF DYNAMICS 

A. System parameters 

In the following we consider Li^ atoms with mass 
M w 1.0 • 10~^^ kg and the scattering length tuned 
down to fls « 2.4 • 10~^°m by using a magnetic Fesh- 
bach resonance; the trap frequencies read w = 27r • 60 
Hz and uj± = 27r • 3.6 kHz. Returning to dimensionless 
units, we obtain the weak coupling g w 0.055. The initial 
spin-polarized equilibrium is characterized by the parti- 
cle number N = Ep = 100 and by the temperature T 
taking values of either T/Ep = 0.2, 1, or 5 correspond- 
ing to the degenerate, intermediate and non-degenerate 
regime. According to these values one finds: the chemical 
potentials ^/Ep sa 1.0, 0.54, -7.2; cloud extensions A « 
28, 40, 89; maximum densities no(0) « 4.4, 3.3, 1.8; and 
maximum mean-field potential strengths no(0)(7 w 0.24, 
0.18, 0.097 that are small compared to the trap frequency, 
being 1 in our units, and extremely small with respect to 
typical single-particle energies max(T, /i) > 100. In ad- 
dition to the temperature, we also vary the spiral wave 
length As and consider either A/ As = 1, 2 or 5 wind- 
ings of the spin spiral within the atom cloud. For these 
conditions, we integrate the time evolution using a Mac- 
Cormack method [l3| and trust results that coincide for 
grid-sizes 300^ and 600^ for a phase space region of linear 
extension sa 2A. 



B. Observation of spin segregation 

On a short time scale ~ 1, the system evolves mainly 
as determined by the harmonic trapping potential. That 
is, neglecting interaction completely for the moment, 
{wq{x,p),w{x^P)^ simply rotates in phase space at con- 
stant angular velocity 1 (in units of the trap) ; each point 




1 t/27i 2 

FIG. 1: (Color online) (a) Wigner functions Wx{^,p) and 
WziXjp) at five times [i-v, as indicated in (b)] during the first 
half cycle; T/Ef = 1, As/ A = 0.5, both x and p range from 
-40 to 40. The motion in phase space is governed by an overall 
rotation at the trap frequency, during which the z-component 
of the Wigner function slowly builds up two domains, (b) 
Time evolution of the averaged spatial polarizations dz and 
(Jxy for T / Ef = 0.2, 1, 5 (thick solid, dashed, and thin solid 
lines, respectively) and for As/A = 1, 0.5, 0.2 (upper triple 
of red, central triple of black, and lower triple of blue lines, 
respectively). While for the shorter wave lengths As /A the 
spatial a;t/-polarization shows rapid collapses followed by pe- 
riodic revivals, the emerging spatial ^-polarization undergoes 
smooth oscillations as a signature of the formation of only two 
domains in phase space. 

of the Wigner function follows a classical circular orbit. 
This behavior can be observed clearly in the first row 
of Fig. [Ija) showing the evolution of during half a 
cycle. This simple dynamics in phase space translates 
into a more involved evolution of the spatial polariza- 
tion n(x), obtained by projecting w{x,p) onto the x- 
axis. n{x) = Jdpw{x,p). Fig. [T](b) shows the averaged 
absolute spatial xy-polarization, 

a^y^j^ Jdx[nl{x)+nl{x)f/\ (19) 

during the first two cycles. A rapid collapse of a^y fol- 
lowed by periodic revivals can be observed, the more pro- 
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FIG. 2: (Color online) Evolution on longer times, for T/Ep = 
1, and different wave lenghts As/A. The l.h.s. shows the av- 
eraged spatial absolute z-polarization, a^, versus time. On 
the r.h.s. the Wigner functions Wx and w, are plotted for five 
instants of times like in Fig. [Ha), but during the 21st half 
cycle (arbitry color scale). Initially increases linearly in 
time with a rate that is controlled by the wave length As/ A. 
Irrespective of the number of windings A/ As (directly visible 
in Wx), Wz developeds two oppositely polarized domains. 

nounccd the larger the number of spiral windings A/Xg 
in the atom cloud. 

During a single cyele of oseillation in the trap the weak 
interaction causes only small deviations from a simple ro- 
tation in phase space. The small modification of the trap 
frequency and the slight anharmonicities caused by the 
scalar part of the mean-field potential Vn-^{{x) are hardly 
noticeable. However, the effect of interaction becomes 
apparent in Wz , being zero initially [Fig. [IJa) , note the 
different color scales for Wx and Wz]- Though still small, 
Wz develops a characteristic pattern induced by the mag- 
netic mean field B^f. Namely, in Wz two domains of 
opposite z-polarization form. This spin segregation in 
phase space corresponds to phase-opposed dipole oscil- 
lations of the t- and the ^-polarized domain in the trap 
(see also Fig. [3]). Remarkably, the spin segregation does 
not reproduce the structure of the initially created spin 
spiral of wave length A^ . The formation of two spin do- 
mains (and only two) is a very robust effect; we find it 
for all pitches of the spin spiral considered here. The up- 
per panel of Fig[TJb) shows the averaged absolute spatial 
z-polarization 

a.E^l jAx[nl{x)f'\ (20) 

As a consequence of spin segregation in phase space, <Tz 
oscillates in time, however, because only two domains are 
formed, it does not feature sharp collapses like (Jxy does 
for small As /A. 

In Fig [2] we present data for longer times, for the orig- 
inal temperature T / Ep = 1 and for three different spiral 
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FIG. 3: (Color online) Spatial densities n^{x) and ni{x) as 
well as z-polarization nz{x) = ^\n^{x) — ni{x)\ at five instants 
of time during the 21st half cycle, for T/Ef = 1 and As/A = 
0.2. Times and parameters correspond to the middle row 
of the r.h.s. of Fig. [21 Spin t and \. particles segregate and 
perform phase-opposed dipole oscillations in the trap. 



wave lengths As/A: From cycle to cycle the spin seg- 
regation becomes more and more pronounced as visible 
from Uz and from the snapshots on the r.h.s. showing Wz 
during the 21st half cycle. The rotation in phase space 
of the two oppositely polarized domains corresponds to 
phase-opposed dipole oscillation of f and \. spins in the 
trap. This behavior is visible in Fig. [3] showing the spa- 
tial densities n^{x) and n],{x) as well as the spatial polar- 
ization nz{x) = ■^[n^{x) — ni{x)] during the 21st cycle. 
While a z-polarization builds up, the spiral spin struc- 
ture in the xy-plane decreases but stays intact (cf. 
during the 21st half cycle shown in Fig. [5]). The total 
density Wigner function wq hardly changes during the 
time-evolution also for longer times (not shown). The 
spin segregation can directly be controlled by the num- 
ber of windings A/ As of the spin spiral within the cloud, 
the more windings the slower the segregation builds up. 
The fastest segregation is observed for A/ As = 1, here al- 
ready after 10 half cycles deviations from a linear increase 
is found and a more complex dynamics sets in (Fig. [5]). 



V. EXPLANATION BY DYNAMICALLY 
INDUCED LONG-RANGE INTERACTION 

A. Effective Description 

In order to give an intuitive explanation for the spin 
segregation, let us describe the system in the rotat- 
ing phase-space frame with the new coordinates x' = 
xcos(t) — psin(t) and p' = pcos(t) -I- a;sin(t) describing 
classical orbits in the trap. In that frame w' {x' , p' , t) = 
w{x{x' ,p' ,t),p(x' ,p' ,t),t) is stationary for vanishing in- 
teraction. However, interaction, represented by the 
mean-field potentials (|18p . is now time-dependent, since 



it is obtaned by projecting onto the x-axis that rotates 
with respect to the new frame. For example, the mag- 
netic mean field B'^-^f{x' ,p' ,t) = Bjni{x{x' ,p' ,t),t) reads 

/oo 
ds w' {x'—s sin(t) , p'+s cos(f), t) . 
-oo 

The time-dependence of the mean field -B[,jf originates 
on the one hand from the rotation of the integration axis 
at trap frequency and on the other hand from the time- 
dependence of the Wigner function w' {x' , p' , t) . In the 
rotating frame, the latter is solely governed by the weak 
interaction and slow compared to the oscillation in the 
trap. We can use this difference in time scales to separate 
the dynamics on short times from that on longer times. 
We assume that a single oscillation in the trap is not 
affected by the weak interaction. This allows us, in turn, 
to integrate out the rapid oscillations in the trap when 
studying the dynamics on longer times where interaction 
does play a role; we approximate 



BnAx',p',t) « B^ff{x,p',t) 
_ -2,9 '•2- 



2tt 



givmg 



dp 



dsw'(^x' — ssm{T),p' + scos(r),t) 



H-r- , ^' — m'jx' 



p2 



{x' + x,p' + p,t). 

(23) 

By averaging over a cycle, we have obtained an 
effective mean-field potential that corresponds to a 
time- independent isotropically long-ranged interaction in 
phase space. By oscillating in the trap, the system dy- 
namically acquires a spatially long-range interaction. 



B. Zero-order semiclassical mean-field interaction 

We can simplify the description further, again argu- 
ing that interaction is weak and the mean-field potential 
small compared to the trap. For the mean-field contribu- 
tion V,^,^^ of the potential t4'm(x) = V^]^^{x) + ^x^ , ap- 
pearing in the infinite series on the r.h.s. of the equation 
of motion (jlip , we truncate the series already after a = 
instead of a = 1 . In the set of equations p5|) , this approx- 
imation corresponds to neglecting the mean-field-induced 
acceleration by dropping terms involving the gradients 
dxBynf and dxVmf. We keep, however, the spin-rotating 
term B^{ x w stemming from the order a = 0. Together 
with cycle averaging, in the rotating phase-space frame, 
we arrive at the effective equations of motion 



w'o{x',p\t) 

w'{x',p\t) ^ B-3{x',p\t)xw'{x',p\t). 



(24) 



The second equation describes the time evolution of the 
polarization field w' in the x'p'-plane. At each point the 
polarization rotates in the mean field B^f such that \w'\ 
stays constant. 






XslA = 0.2 X,IA = 0.5 X,IA=l 



FIG. 4: (Color online) Cycle-averaged rate of the creation 
of z-polarization computed for the initial state, [i3J5^(0) x 
iy'(0)]z, for T/Ef — 1 and different spiral wave lenghts As/A 
(arbitrary color scale). It explains the formation of two op- 
positely polarized domains in Wz that is visible in Fig. (2] 



C. Growth of 2-polarization 

The spin segregation observed numerically can now 
be explained by first-order time-dependent perturbation 
theory, predicting according to Eqs. (|24p initially a linear 
growth of the z-polarization, 

w,{x',p',t) ~ [B,ff(x',p',0) X w\x',p\0)]A (25) 

as we can observe it on the l.h.s. of Fig. [2] Deviations 
(22) from the linear growth appear as soon as w'^ be- 
comes comparable to as visible in the lower left plot 
of Fig. [3 Figure H shows the rate [-B,'ff(0) x w'{0)]^ com- 
puted for the intermediate temperature T / Ep = 1 and 
different spiral wave lengths A^/A. Notably [BJ^(O) x 
w'{Q)\z always shows a pattern with two oppositely po- 
larized domains, irrespective of the number of windings 
A/ As [iBl- This explains the previously observed seg- 
regation of t and ^ polarization. The formation of two 
domains only can be understood as follows: 

According to Eq. the spin polarization w' at 

a given point [x' ,p') in phase space feels a magnetic 
mean field that mainly depends on the polarization found 
in phase-space areas close by. Within the phase-space 
neighborhood of {x' ,p'), in turn, phase-space areas show- 
ing the largest polarization perpendicular to w'{x',p') 
contribute most. Since for the initial state \w'\ ~ Wq/2 
increases towards the origin x' — p' — 0, at a given point 
{x',p') the effective magnetic mean field B^ff is domi- 
nated by the polarization found when slightly moving 
along the direction of the spiral towards the origin. On 
one side of the spiral this results always in the creation 
of a positive z-polarization, on the other side always to 
the creation of a negative z-polarization. This explains 
the creation of two domains. Moreover, the strength of 
the local mean field depends on the spatial variation of 
the spin density \w'\ compared to the spiral wave length 
A^.; the smaller Ag/A the slower the spin segregation, as 
observed numerically in Fig. 

We have identified the mechanism underlying the ob- 
served spin segregation. Obviously, the phenomenon 
does not depend on the sign of the spin-spin coupling. 
We have checked numerically that it is equally observ- 
able for attractive interaction, giving reversed polar- 
izations. Moreover, it can also be expected for non- 
condcnscd "spin-1/2" bosons which are equally described 
by Eqs. p5|) . but with the exchange interaction giving 
rise to a spin-coupling of opposite sign. 
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VI. EXPERIMENTAL SIGNATURES 

One can measure the dynamical spin segregation by 
state-sensitive absorption imaging either in situ (as in 
the experiment by Du et al. [1]) or after a time of flight. 
In the latter case one can also use Stern- Gerlach sep- 
aration to distinguish f and ], particles. An in situ 
measurement gives the spatial distributions nm{x) = 
Jdpwmmix,p) of both Spin states m the images 

after a time of flight reveal their momentum distribu- 
tion fimip) = J<ixWmmix,p). One can then determine 
the z-polarization in space, nz{x) = ^[n-f-(a;) — n^(x)], 

and momentum, nz{p) = ^["■tC?') ~ '^i(p)]- shown 
in Fig. [3l the dynamical spin segregation corresponds 
to phase-opposed dipole oscillations of f and J, spins in 
space. The momentum distributions will show the very 
same behavior, but shifted in time by the quarter of a 
cycle (because momentum densities are obtained by pro- 
jecting the Wigner function onto the p-axis) . 

VII. CONCLUSIONS 

The phenomenon of dynamical spin segregation pre- 
dicted here is different from the effect observed by Du et 
al. described in Refs. [I,!!]. In their case no spin spiral is 
created initially, instead an inhomogeneous external mag- 
netic field is present throughout, leading eventually to a 
spherical symmetric spin segregation in phase space be- 
tween an inner core and an outer shell. The phenomenon 
described here also differs from the physics of the spin- 
wave instability investigated by Conduit and Altman Q . 
They consider the same spiral spin structure as initial 
state, but strong repulsive interaction. In contrast to the 
dynamical spin segregation into two counter oscillating 
domains found here, crucially depending on the presence 
of the trap, their spin-wave instability leads to spatial 
(non oscillatory) domain formation, not requiring a trap 
and with the domain size controlled by the spiral wave 
length. 

The system's dynamics described here can be called 
self driven. The transformation to the co-rotating 



frame in phase space corresponds to the transforma- 
tion to the Dirac picture on the full quantum many- 
body level, where the quadratic Hamiltonian i^o = 
H — Hint constitutes the unperturbed problem. In the 
Dirac picture, the time evolution is generated solely by 
the time-dependent interaction Hamiltonian H^^^^{t) ~ 

exp(— iifoi/ft)-ffint exp(iffoV^)- Thanks to the equidis- 
tant ladder spectrum of the harmonic trap, it is time 
periodic, H^^^^{t + T) = -fffnt(^) with T = 27r/w, like the 
Hamiltonian of a driven system. This additional sym- 
metry has strong consequences for the dynamics. It al- 
lows us to find a time-independent effective description 
iJgff = ^ /o^dt TJj'Jjj (i) for the dynamics on longer time 
scales that does not depend on the details of the short- 
time dynamics. For the system described here, the ef- 
fective description, in form of the mean field potential 
(|23p , contains a spatially long-range interaction (isotropic 
in phase space) that the original Hamiltonian did not 
possess and that explains the spin segregation observed. 
The system has dynamically acquired novel properties. 
A similar situation is found for example for interacting 
particles in tilted lattice systems with the single-particle 
spectrum given by the Wannier-Stark ladder [ij. The 
separation of time scales found here resembles also the 
physics of driven many-body systems as it has been stud- 
ied in lattice systems subjected to off-resonant exter- 
nal driving. These systems are equally described by an 
approximate effective time-independent Hamiltonian on 
long times [l5| . 

The dynamical spin segregation is ~- like itinerant ferro- 
magnetism - caused by exchange interaction. However, 
while the Stoner transition to a ferromagnetic phase is 
an equilibrium effect requiring fairly strong interaction 
(as well as spatial dimensionalities larger than one), the 
robust effect described here happens far from equilibrium 
and does not need strong interaction but, instead, suffi- 
ciently long times to build up. 
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